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ATTITUDE REORIENTATION OF SPACECRAFT BY MEANS 


OF IMPULSE CONING 

C. William Martz 
Langley Research Center 


SUMMARY 

An investigation has been conducted to determine minimum "maneuver costs" 
for attitude reorientation of spacecraft of all possible inertial distribution 
over a wide range of maneuver angles by use of the impulse coning method of 
reorientation. Maneuver cost is proportional to the product of fuel consumed 
and time expended during a maneuver. Assumptions included impulsive external 
control torques, rigid-body spacecraft, rest-to-rest maneuvers, and no distur- 
bance torques. Also, coning maneuvers were constrained to have equal initial 
and final cone angles. 

Maneuver costs are presented for general reorientations as well as for 
spin-axis reorientations where final attitude about the spin axis is arbitrary. 


INTRODUCTION 

Impulse coning is a method of attitude reorientation in which impulsive 
torques are used to initiate and (later) to terminate free precessional motion 
of a body. Between torque impulses, the spinning body coasts undisturbed 
through some desired attitude change. Impulse coning has potential application 
to all maneuvering spacecraft equipped to generate directed external torques. 

In past years, attitude reorientation by means of impulse coning has been 
examined by several authors (refs. 1 to 4). Reference 1 suggests that large- 
angle reorientations of axisymmetric bodies be accomplished with a series of 
"small-angle" precessional motions. Reference 2 presents a two-impulse coning 
scheme for large-angle spin-axis reorientations of axisymmetric bodies. Spin- 
axis or pointing-axis reorientations can be defined as general reorientations 
with final attitude about the spin or pointing axis arbitrary. Reference 3 
extends the spin-axis reorientations of reference 2 to asymmetric bodies. How- 
ever, precession angle is restricted to 180° and spin-axis inertia must be 
larger than inertias about transverse axes. Reference 4 introduces a statis- 
tical "average cost" for general reorientations over a specified maneuver-angle 
range and favorably compares two-impulse coning costs with costs of other 
reorientation methods. Results are limited to elongated axisymmetric bodies, 
however . 

The present paper extends previous results to include two-impulse coning 
costs for general reorientations of axisymmetric foreshortened bodies. Also, 
two-impulse coning costs are obtained for spin-axis reorientations of asymmet- 
ric bodies of all possible inertial distribution. Finally, maneuver costs are 



determined for general reorientations of asymmetric bodies. In this category, 
maneuvers consist of a coning motion followed by an impulsive spin maneuver to 
obtain the desired spin attitude. 

In the present investigation, reaction thrusters are assumed to initiate 
and terminate all coning motions with impulsive torques of negligible duration. 
Other assumptions include rigid-body spacecraft, lack of disturbance torques, 
and rest-to-rest maneuvers (zero initial and final attitude rates). 

The product of total maneuver impulse (proportional to fuel consumed) and 
total maneuver time is determined for the various solution paths available for 
each maneuver. This product is then nondimensionalized by spacecraft inertia 
about the intermediate principal axis and is referred to as the cost function 
for that maneuver. Since total impulse and maneuver time are inversely propor- 
tional for impulsive maneuvers, results are independent of each of these quan- 
tities. The desired or optimal solution is defined as the solution associated 
with the smallest cost function for each maneuver. Since it was assumed that 
all reorientations within a given range of maneuver angle are equally probable, 
optimal cost functions for a large number of statistically representative 
reorientations within a given maneuver range were averaged to represent reori- 
entation costs for that maneuver range. Maneuver ranges of 0.00873, ir/4, tt/2, 
and 17 radians were used. 


SYMBOLS 


F 

FLAG 

H 


> ^y » ^z 
Hxy 


Ixj^yj^z 


m 

n 


nondimensional asymmetry factor defined by equation (17) 
indicator of nutation class (see eq. (12)) 

total angular momentum, defined as positive for ij; > 0 and negative 
for ij) < 0 

angular momentum components along x-,y-,z-axes, respectively 

angular momentum in transverse plane, + Hy^ 

transverse inertia for roll-symmetric bodies 

principal body inertias about x-,y-,z-axes, respectively; ly is 
intermediate moment of inertia 

mass element 

= 0 , 1 , 2 , 3 , . . . 


Qp product of impulse and time for thruster prepositioning 

Qg product of total impulse and total time divided by ly for spin 

maneuver (see eq. (38)) 
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Q2 

Q2,A 

Q2,T 

^2 , swivel 
Q3 

Rl »R2>^3 


^2 , max 
sgn ( ) 

T 

t 

tf 

ts 

X,Y,Z 


x,y,z 


Xd 


xf 

f 

Xq 

y’ 

A(|) 


e 


0m 

ip,e,4) 


product of total impulse and total time divided by ly for coning 
maneuver (fixed thrusters) 

axial component of Q2 (see eqs. 36)) 

transverse component of Q2 (see eqs. (36)) 

product of total impulse and total time divided by ly for 
coning maneuver (thrusters on spin axis swiveled) 

product of total impulse and total time divided by ly for three- 
impulse coning- spin maneuver 

sequential Euler angle rotations about z-,y-,z-axes, respectively, 
relating body axes before and after reorientation (see fig. 4 ), 
R2 also referred to as maneuver angle 

maneuver range, radians; 0 ^ R2 = R2,max 

sign of quantity 

kinetic energy (see eq. (23)) 

time, sec 

coning maneuver time, sec 
spin maneuver time, sec 
inertial coordinate axis system 

body principal coordinate axis system; unless otherwise specified, 
z is spin axis 

X(j-axis translated from mass center along Zf- axis 
Xf-axis translated from mass center along Zf-axis 
Xg-axis translated from mass center along Zg-axis 

y-axis translated from mass center along Zg-axis after R-| rotation 
change in spin angle during coning, (t>f - (j>g 

error in spin angle at termination of coning maneuver (see fig. 4 ) 

extremum value of 9 ( see text following eq . ( 19 )) 

rotation sequence of Euler angles about x-,y-,z-axes, respectively, 
relating inertial and body coordinate axis systems 
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Euler rates 


4>g "spherical" values of defined by equation (6) 

a)x,Wy, 0)2 body rates about x-,y-,z-axes, respectively 
Subscripts: 

av average value (with respect to time) during maneuver 

d desired value 

f value at conclusion of coning maneuver 

max maximum value 

min minimum value 

0 value at start of coning maneuver 


ANALYSIS 

Equations of Motion 

Two coordinate axis systems used in the analysis of the impulse coning 
method of attitude reorientation are illustrated in figure 1 . The body prin- 
cipal axis system x,y,z is related to an inertial axis system X,Y,Z by an 
Euler rotation sequence about body axes z,y,z, respectively. Note 

that 4^ and Z are chosen to be in the same direction as the total angular 
momentum H. Angular momentum H is established for each reorientation maneu- 
ver by the initial torque impulse of the thrusters. Reaction thrusters are 
assumed to be fixed along the x,y,z body axes to produce torques about the 
y-,z-,x-axes, respectively. Without loss of generality, 6 is restricted to 
0 < 6 < TT/2. 

From figure 1(a), the relationships between inertial rates and body rates 
Wx> ^y> and are 

(Oy sin 4) - Wjj cos <t> 

i> = T 

sin 0 

0 = Wy cos <t> + Wjj sin 4> / ( 1 ) 

4> = 0)2 - 4* cos 0 
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The angular momenta along the body axes (from fig. 1(b)) are 


H cos 6 = 

H sin 6 sin (|) = 
-H sin 9 cos <|) = 


( 2 ) 


Combining equations (1) and (2) yields the Euler rate equations governing 
the coning motion, which are 


= H 


/sin^ 4* cos^ 4>^ 


/I 1 \ I 

6 = H — - — I sin 6 sin (j> cos <)> / 


cos 0 

<j) = H - 4) cos 0 


(3) 


Axisjnnme^ric Bodies 


The Euler rate equations (eqs. (3)) are easily integrated for axisymmetric 
dx = ly = I) bodies with the result for 4>o = ^ expressed as 


H H 

4^ = — t and 4’f = — tf 
I I 


9 = 0Q = = Constant 


4) = 4>0 + 4^ cos 0n 1 


and 


A4> = 4>f - 4>o = 'l^f Gos 0o( "I 


(4) 
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This solution is the free precessional motion of a spinning body as illustrated 
in figure 2. The specified initial and final locations of the body z (spin) 
axis are elements of a space cone connected by a coning or precession circle. 
The total angular momentum H and the Z inertial axis are coincident with 
the cone axis. 

An impulse coning maneuver is performed as follows: The body, initially 

at rest, is acted upon by an external control torque impulse of negligible 
duration which causes it to spin about its z-axis and process about the Z-axis. 
When the spin axis reaches the desired final location, a second external torque 
impulse terminates the precessional motion. If the body attitude about the 
spin axis 4» is as desired, the second impulse also terminates the spinning 
motion. Otherwise, the second impulse will include a component along the spin 
axis to initiate a corrective spin maneuver which is terminated by a third 
impulse. 

The axial and transverse momentum components generated by the initial 
impulse are labeled o ^xy,o figure 2. Also, the momentum compo- 
nents which terminate the motion are shown as H2^f and The initial 

body-axis system XQ,yg,ZQ is included to illustrate the orientation of H 
in body coordinates. 

It should be mentioned that candidate coning solutions to equations (3) 
must include negative as well as positive precession rates and spin rates. For 
example, the two pairs of alternative paths illustrated in figure 3 are similar 
geometrically but differ in the sign of 4’ and (j) and therefore have unlike 
final spin attitudes. 

The spherical geometry of a generalized reorientation maneuver R-|,R2,R3 
is shown in figure 4. By definition, Ri,R2,R3 is any specified set of 
sequential Euler angle rotations about the z,y,z body axes, respectively, 
which reorient the body from its initial attitude to some final desired 
attitude . 

Parameter relationships for the coning motion can be determined from the 
geometry of figure 4. The oblique spherical triangle with sides 0 q, R 2, 

and 0f is bisected into two right spherical triangles with mirror symmetry. 
From the leftmost triangle of figure 4(a) 


sin Sq 


sin (R2/2) 
sin (4»f/2) 


(5) 


and 


tan (R2/2) 

sin (4's/2) = sgn 

tan Qq 


( 6 ) 
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where sgn ilJ is used to control the sign of ij^g. At the start of coning (see 
fig. 4) 


'I's TT 

<|)o + R-| = -(2 - sgn i|)) 


(7) 


At the end of coning 


(|)f + 



e 



( 8 ) 


Equations (7) and (8) are combined to yield the following equation for the 
error in spin angle: 


e = Rl + R3 - vjig - A(t) + irCegn ip - D 

and because of its cyclic nature, e can be written 

e = Rl + R3 - \j)g - A<t) ± 2Tm (n = 0, 1, 2, . . .) (9) 

where the value of n is adjusted to minimize the absolute value of e to 
eliminate unnecessary maneuver costs. 

Equations (5), ( 6 ), (4), (7), and (9) govern the coning solution and 
relate the variables tfi 9, 'I'si <t>o> ^>^8 e, for a given maneuver 

(Ri,R 2,R3) and inertial distribution (I^ = ly Ij;). Also, these variables and 
constants combine to generate a cost function for each maneuver as reported in 
the section entitled "Reorientation Costs." Because the solution is under- 
defined (more variables than independent equations) it was necessary to perform 
a trial -and-error search (search variable 'l^f) to determine the solution with 
the smallest cost function for each maneuver. All such solutions were found 
to occur at e = 0. That is, two-impulse coning not only is sufficient for all 
reorientations of axisymmetric bodies but also results in minimum reorientation 
costs. 


Asymmetric Bodies 

For asymmetric bodies (ix > ly > Iz), the Euler rate equations (eqs. (3)) 
are nonlinear, and analytical solutions are not available except in the form 
of elliptic functions. A more useful approach is the application of Poinsot's 
construction which is presented in reference 5. This method qualitatively 
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analyzes the free motions of rigid bodies and indicates that the motion of asym- 
metric bodies is similar in many ways to that of symmetric bodies. A notable 
exception is the nutational or nodding motions of the z-axis as it precesses 
about the momentum vector. Thus, for asymmetric bodies, 0 is oscillatory as 
are the time histories of the Euler rates. 

In order to exploit the energy-momentum relationships of the Poinsot method 
and determine quantitative information regarding the free motion of asymmetric 
bodies, it was decided to restrict solutions to the class for which 6f = 0 q. 
This restriction not only furnished an additional equation in 0, but imposed 
a symmetry of motion on all solutions which in turn constrained (j) and elimi- 
nated the need for trial-and-error solutions. 

Figure 2 illustrates nutational motion during the coning maneuver by the 
curves labeled "Inner path" and "Outer path." That is, there are two separate 
solutions for each maneuver which satisfy all initial conditions as well as the 
imposed terminal condition 0f = 0 q. These solutions are shown quantitatively 
in figure 5 which presents normalized 0-4) trajectories obtained by numerically 
integrating equations (3) for an example maneuver. Because of the terminal 
condition (0f = 0 q)» all such trajectories are characterized by a symmetry of 
motion about the midtrajectory point. It can be shown with equations (3) that 
spinning motion associated with inward nutation is antisymmetric about the mid- 
trajectory point <1) = 0 ± TTn; this leads to the following roll-constraint equa- 
tion for inward nutation: 


<t>f + 4 )q = 0 ± 2iTn 


(n = 0, 1, 2, . . .) (10) 


The last term of the two equations in the preceding sentence derives from the 
symmetry of the inertia ellipsoid (ref. 5) about the body spin axis and the 
periodic nature of <|). If 4>o is changed by increments of ±2iT, ±4 it, ±6tt, etc., 
before the initial coning impulse, the coning motion will be unchanged. Also, 

4>o can be changed by increments of ±tt, ±3Tf, ±5 it, etc., without changing the 
coning motion providing the transverse impulse torques are applied in the 
opposite direction. These alternatives are implied by the last term of 
equation (10). 

Referring again to figure 5, the path labeled "Outward nutation" is part 
of a trajectory set for which the spinning motion is antisymmetrically centered 
about 4> = TT/2 ± TVn causing the roll-constraint equation for outward nutation 
to be 


4>f + <t>Q = ir ± 2iTn 


(n = 0, 1, 2, . . .) (11) 


The last term is included for reasons just mentioned in connection with equa- 
tion (10). For manipulative convenience, equations (10) and (11) are combined 
into a single constraint equation which is 


8 


(l>f. + (J>Q = tt(FLAG) ± 2m 


(n = 0, 1, 2, . . .) 


( 12 ) 


where FLAG = 0 for inward nutation and FLAG = 1 for outward nutation. This 
roll constraint precludes two-impulse coning reorientations of asymmetric 
bodies for the condition 6f = 9 q except for a few maneuvers which satisfy the 
following combination of equations (7), (8), and (12) for e = 0 ± 27®!; 


R3 - Rl = ir(FLAG - 1 ± 2n) (n = 0 , 1 , 2 , . . . ) (13) 


All other coning reorientations of asymmetric bodies require a follow-up spin 
correction maneuver. 

Another relationship pertinent to the coning motions of asymmetric bodies 
concerns precession time and is written 


^f - 4^0 

tf = — ^ 

4^av 


(14) 


where = 0 and (jjav is determined by integrating the first of equations (3) 
to obtain 


or 




Ix + ly 

21xly 


ly - Ix sin (2<|)f) - sin (2(1 )q) 
21xly 2(A(1)) 


(15) 


Equation (15) represents the average ^ with respect to <{) over the coning 
maneuver and thus an approximation to the average 4) with respect to time 
required by equation (14). However, the error of this approximation was found 
to be negligible for the maneuvers of the present paper. For example, this 
error is less than one part in ten thousand for the trajectories of figure 5. 
Combining equations (12), (14), and (15) and simplifying yield 


l|)flx 

hTf) 


( 16 ) 
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where 


ly p sin (A(|)) 


[l - 2 (FLAG)] 


An independent source for A(|) is obtained by integrating the third of equa- 
tions (3) to yield 


ptf ^ ptf/fl 

A(l)zJ 4)dt=J \^- 


cos 6 . \ 

- ^ cos 0 jdt 


Replacing cos 0 with its average value during the coning maneuver and inte- 
grating give 


HtfCcos 0)av 

A(|) ipf(cos 0)ay 

It; 


The last term is an approximation to the original integral and is not accurate 
for maneuvers having large variations in 0. However, 0 variations as moni- 
tored for optimum reorientations were small without exception and, by using 
the following approximation for (cos 6)av» 




cos 0 d0 2 


cos 0 d0 


(cos 0). 


0f - 0o 


2(0„, - 0r.) 


sin 0m - sin 0 q 
6m - 00 


equation (18) has been found to compare closely with numerically integrated 
results. In equation (19), 0m is the limiting or extremum value of 0 as 
given by equation (31) or (32) or (33) or (34), depending upon the polhode axis 
and whether the nutation is inward or outward, as explained in the section 
entitled "Nutational Motion Limits.” 

Combining equations (16), (18), and (19) gives 


A(|) = \l;f 


sin 0m - sin 0o\/Ix/Iz 

0m - 6o A F 
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Equations (12), (7), (6), (31) or (32), (5), (20), and (9), which govern 
variables <l>o. 'Ps, Qq* r’espectively , were solved 

iteratively in the order specified by the Newton-Raphson method for each given 
maneuver and inertial distribution. The amplitude of A<]) used to start the 
iterative solution was set initially at 0.000001 and successively increased 
to 0.5, 1, and 2 until convergence occurred. Four separate solution paths 
(one for each of the 6jjj values referred to in connection with eq. (19)) were 
available for both positive and negative precession rates - a total of eight 
possible solutions for each reorientation. All such solutions were determined 
and the one associated with minimum maneuver costs was retained for computing 
the "average cost" of reorientations presented for each maneuver range. The 
listing and documentation of computer program ASYMCON, used to determine these 
solutions, is presented in appendix A. 


Polhode Axes 

Poinsot's geometric solution, described in references 5 and 6, equates 
the free rotational motion of a rigid body to the motion of the "ellipsoid of 
inertia" for that body as it rolls on the "invariable plane." The point of 
contact between these two surfaces traces out a "polhode" curve on the ellip- 
soid of inertia. The axis surrounded by the polhode curve, or "polhode axis," 
is the body axis which processes about the momentum vector during a coning 
maneuver. The polhode axis, which can be an axis of minimum moment of inertia 
or an axis of maximum moment of inertia but not an axis of intermediate moment 
of inertia, is determined by the initial motion and inertial distribution of 
the body. 

It is necessary to determine the polhode axis for each coning maneuver 
before the angular limits of the nutational motion can be determined. Consider 
first the class of bodies for which 


Ix > ly > Iz 


Reference 6 indicates that the condition for polhode about the z-axis is 


h2 

ly > — > Iz 

y 2T 

where or (Momentum) is expressed 

h 2 = dx^x)^ + + dz^^z)^ 


( 21 ) 


( 22 ) 


11 


and 2T, or 2 (Kinetic energy), is expressed 




( 23 ) 


Equations (22) and (23) > with body rates substituted from equations (2), 
are combined with that part of inequality (21) concerned with intermediate 
inertia ly, and the condition becomes, for polhode about the z-axis, 


Ixdy - Iz^ 

tan2 0 cos2 0 < (I^ > l„ > l„) 

Izdx - ly) X y z 


where 9 and <t> are chosen to be the initial conditions 0 q and 4>o* The 
condition for polhode about the x-axis is 


h2 

Ix > — > ly 
2T ^ 


(24) 


Again the momentum and energy expressions, with body rates from equations (2), 
are substituted into the part of inequality (24) involving ly, and the condi- 
tion becomes, for polhode about the x-axis. 


Tx(Ty “ I 2 ) 

tan^ 0 cos^ <|) > - — dv > Iv ^ Iz) 

Izdx - ly) X y z 


Now consider the remaining asymmetric bodies which are in the class 


Ix < ly < Iz 

The condition for polhode about the z-axis is 
r2 

I^> —> 


Substituting the momentum and energy expressions and body rates into the ly 
part of inequality (25) results in the condition 


9 o 5 ^xdy - Iz) 

tan2 0 cos2 ^ < (I„ < I„ < I„) 

Izdx - ly) X y z 
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for polhode about the z-axis. Finally, the condition for polhode about the 
x-axis is 


h2 

Iv > — > Ix (26) 

y 2T 

Again, substituting the momentum and energy expressions and body rates into the 
part of polhode condition (26) involving ly yields the condition for polhode 
about the x-axis 


Ixdy - Iz) 

tan2 0 cos2 (j) > (Iy < Iv < I-) 

Ifl-I) xyz 


In summary, polhode is about the x-axis if 


Ixdy - Iz) 

tan2 0 cos2 4) > (27) 

^z(^x ~ ^y) 


Conversely, polhode is about the z-axis if 


Ixdy - Iz) 

tan2 0 cos2 <|) < (28) 

^z(^x “ %) 


Nutational Motion Limits 


As the spin axis of a rigid asymmetric body precesses about the angular 
momentum vector in free rotational motion, a nutation or nodding of the axis 
occurs which causes the half-cone angle 8 to oscillate between some maximum 
and minimum values. These limiting values are derived in reference 5 and pre- 
sented along with an analysis of the motion. The limits are needed in the 
present investigation to determine the mean value of the half-cone angle during 
the coning maneuver for use with equation (19). 

Reference 5 shows that for a polhode about the z-axis and either 
Ix ^ ly > Iz or Ix ^ ly < Iz 


tan^ ®inax 


Iy(H2 - 2TIx) 

Ix(2TIy - h2) 


(29) 
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and 


tan2 


Ix(h 2 - 2TIz) 

Iz(2TIx - h2) 


(30) 


where Oniax and ©niin are the limiting values of 6 over a full precession 
cycle. The coning maneuvers of the present paper use only a small part of this 
cycle, however. Thus, only Omax ia needed for outward nutations and only 
®min ia needed for inward nutations. 

These motion limits can be expressed as a function of inertial distribu- 
tion and Euler angles by substitutions for the momentum and energy (eqs. (22) 
and (23)) and the body rate (eqs. (2)). The following equations result for 
polhode about the z-axis: 


a 

°max 

and 

®min 


tan“^ 



lydz - Ix) + Ixdz - ly) tan2 ^ 

j 

l 2 dx - ly) + Ixdz - Iy)/(tan2 6 cos2 


tan“^ 



lydz - Ix) + Ixdz - ly) tan2 (j> 

y 

Izdy “ Ix) tan2 <t> + lydz " Ix)/(tan2 6 cos2 (|)) 


(31) 


(32) 


where <t> and 9 are chosen to be evaluated at the instant coning starts. 

Thus far, limits have been determined only for motions characterized by 
polhodes about the z-axis. Consider now the alternative maneuvers with 
polhodes about the x-axis. The motion-limit formulas of reference 5 can be 
related to Qmax ®min only when 9 is measured from the polhode axis. 

Rather than rederive the equations of motion based on 9 being measured from 
the x-axis, a two-step alternative was employed to use the existing equations 
of motion for maneuvers having polhodes about the x-axis. The first step is 
to interchange the values of and The second step is to replace the 

original Euler rotations about the z-, y-, and z-axes with a newly generated set 
of Euler rotations about the -x-, y-, and -x-axes designed to produce the iden- 
tical body reorientation. The new set of Euler angles was determined by the 
method of reference 7 and implemented in program ASYMCON of appendix A. With 
these changes, the coning geometry is transferred, in effect, from polhode 
about the z-axis to polhode about the x-axis. 

Reference 5 indicates that for a polhode about the x-axis and for either 
Ix < ^y = ^z- ^x ^ ^y = ^z> limiting values of 9 can be determined 

from equations (29) and (30) with substitutions of equations (22), (23), 
and (2), providing Ix and ly are interchanged. The resulting motion limits 
for polhode about the x-axis are 
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( 33 ) 


®max = 


-1 




IjjCly - Ix) fcan^ (}) + Ivdz - ly) 


y'-^z 


Izdy - Ix)/(tan2 0 cos^ - l^) 


and 


®min ” ^ 


i 


Izdy - ly) tan^ (j) + lydz ~ ^ 


lydz - Ix)/(tan^ 0 cos^ <j>) + lydz “ ^y^ tan^ (j) 


( 34 ) 


where 0 and <1) again are chosen to be evaluated at the instant coning starts. 
Note that these motion limits for polhode about the x-axis can be obtained from 
equations (31) and (32) (motion limits for polhode about the z-axis) simply by 
interchanging ly and Iz- 


REORIENTATION COSTS 

As mentioned in the Introduction, a cost function equal to the product of 
total maneuver impulse (proportional to fuel consumed) and total maneuver time, 
nondimensionalized by spacecraft inertia about the intermediate principal axis, 
is determined for the various solution paths available for each maneuver. The 
desired or optimum path is defined as the path associated with the smallest 
cost function for each coning maneuver. Since total impulse and maneuver time 
are inversely proportional for impulsive maneuvers, results are independent of 
each of these quantities. 

It is assumed that all reorientations within a given maneuver range R 2 ,max 
are equally probable. Therefore, reorientation costs were sampled within each 
given maneuver range by systematically computing optimum cost functions for 
combinations of Ri, R2, and R3 at equally spaced intervals throughout the 
sample space limits -tt ^ R-1 < ir, 0 ^ R2 S R2 max> -it ^ R3 S it where 

R2 max " 0.0087, ■JT/ 4 , iT/ 2 , and it. These optimum cost functions were then 
weighed (according to likelihood of occurrence) and averaged as follows: 


Average cost = 


Z( Optimum cost x sin R2) 
I (sin R2) 


In addition, the standard deviation of the optimum cost functions was deter- 
mined with the relationship 


Standard deviation = 


Ie |j( Optimum cost) 2 x sin R2] 


E(sin R2) 


(Average cost)^ 
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Cost Equations 


The impulse (and thus momentum) required to initiate and terminate coning 
maneuvers has two components (see fig. 2); the transverse component 


^xy,o = ^xy,f “ ^ ®o 


and the axial component 


Hz,o = = H cos 0Q 


For reaction thrusters fixed along the x, y, and z body axes to produce 
torques about the y-, z-, and x-axes, respectively, the cost function 
(Total impulse x Total time)/ly for the two- impulse coning maneuver is 


Q2 = Q2,t + Q2,A 


(35) 


where the transverse and axial components are 
1h| sin (6o)tf 

Q2,t = ; (I sin + [cos 4>ol + I sin ^f\ 

2|h| cos (0Q)tf 
«2,A = 

■^y 







(36) 


If the thruster pair fixed on the spin axis to produce torques about the x-axis 
is allowed to swivel about the spin axis, the components of Q 2 x °sn be mini- 
mized (reduced about 21 percent on the average) at the expense of a relatively 
small thruster prepositioning term Qp with the result expressed as 


|2H(sin Qo^if 2 H(cos ©o^ifl 
Q 2 , swivel = ; + Qf 


Because the thruster prepositioning term introduces variables which make cost 
comparisons between swiveled and fixed thrusters somewhat arbitrary, maneuver 
costs for swiveled thrusters will not be pursued beyond this equation. 
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For axisymmetric bodies (Ix = ly = I), the precession time tf is deter- 
mined from the equation 


tf 


ipf - 4)0 % - 0 



H/I 


and the cost function becomes 


Q 2 = l4^fl sin Sp^lsin 4>o I + loos <l>o I + lain <|)f| + |cos + 2|il'fl cos 6© 

(37) 

Turning now to asymmetric (ix > ly > Iz) todies, it is necessary because 
of the assumed terminal constraint 9f = Qq to combine the coning maneuver 
with a final spin maneuver to attain the desired roll attitude. The resulting 
three-impulse cost function is 


Total impulse x Total time f^2 Qs\ 

T 


where the cost function for the spin maneuver is 


2 e I, 


Q«a = 


- 0.5 1 + 


I4^f^l\ /tg) 


(38) 


The last term of Qg derives from the fact that axial momentum at termination 
of coning can be used to reduce the momentum required for the spin correction 
maneuver provided the momenta are in the same direction. 

From equations (36) and (16), 


Q2,A = 


Q2,t = 


2Ixl4'fl cos 


lyF 


IxliPfl sin 


) (39) 


lyF 


(1 sin (|>ol + I cos <t>ol I sin 4>fl + I cos 4>f I) 
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For minimum Q3 


(IQ3 2 le|i 2 /tf \2 

= 0 = Q 2 — 

d(ts/tf) ly \tg/ 

d2Q3 2leli2/tf\3^ 

d(tg/tf)2 ly \ts/ 


0.5 1 


^ 'I'fE 






J 


(40) 


The time constraint introduced by equations (40) causes the minimum three- 
impulse cost function to be 


QS.min 



(41) 


where for ipfS 0, 


tf |Q2,T + Q2,A 

tg ^ 2I £ I I^/Iy 


and for 'I'fE > 0, 


tf j Q2 ^t 

t^ ■ \|2|€|lz/Iy 

or 

tf Q2,A 

~s ■ 2le| Iz/Iy 


(42) 


(43a) 


(43b) 


whichever is larger. This singularity in tf/tg is caused by the minimum cost 
requirement that the spin correction momentum be larger than the axial momentum 
at the end of coning. 

For application of equation (4l), the iterative solution discussed just 
subsequent to equation (20) was used with equations (38), (39), and (42) 
or (43). 
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Spin-Axis Reorientation Costs 


All spin-axis reorientations were accomplished with two-impulse coning 
maneuvers. For roll-symmetric bodies, the cost function is given by equa- 
tion (37). For asymmetric bodies, the cost function is specified by equa- 
tions (35) and (39). 


RESULTS 

The results of this investigation are primarily the costs of reorientation 
by the impulse coning method. These are presented in the form of average cost 
as functions of 1 - Ig/Iy and R 2 max bodies of all possible inertial 

distribution which are identified in appendix B and presented in figure 6. 
Results for general reorientations are presented first, followed by results for 
spin-axis reorientations and certain cost comparisons. Spin-axis reorienta- 
tions are defined as reorientations with final attitude about the spin axis 
arbitrary. 

Figure 7 shows the effect of inertial distribution on average cost of 
reorientation for maneuver ranges of tt, n/2, 11/4, and 0,00873 radians. The 
lowest curve (I^ = ly) of each plot represents the family of axisymmetric 
bodies. This curve contains minimal two-impulse coning costs and is seen to 
be in agreement with similar' results available from reference 4 for elongated 
bodies with R 2 max = average cost of maneuvering is seen to decrease 

as bodies become less foreshortened or more elongated. 

The remaining curves in figure 7 are for asymmetric bodies and represent 
coning trajectories restricted to the terminal condition 0f = 6 q. With few 
exceptions, these results for asymmetric bodies do not share the lower maneu- 
ver costs of axisymmetric bodies. Comparable reorientation costs for asymmet- 
ric bodies are greater for two reasons. First, the nutational motion prevents 
many potentially low-cost trajectories from meeting the terminal condition 
6f = 0Q. Second, as mentioned in connection with equation (13), almost all 
coning trajectories of asymmetric bodies are unable to satisfy the roll- 
constraint equation without a follow-up spin correction maneuver and increased 
costs . These reasons are applicable to asymmetries as small as one part in a 
million as indicated by a comparison of the = ly and * ly curves. 

In this comparison, the abrupt jump in maneuver costs due to the slightest 
asymmetry is real but, in a practical sense, could be attributed to excessive 
accuracy requirements on the 0 and (j) terminal constraints. In other words, 
if these terminal constraints were relaxed to x parts error per million 
parts, then bodies with asymmetries in I^ or ly up to x parts per million 
could be reoriented by the same impulses used for the Ix = ly results. Thus, 
the Ix = ly results can be used for small asymmetries to the extent that 
errors can be accepted in the terminal values for 0 and <|). For larger asym- 
metries, the other asymmetric curves can be used. Average costs in figure 7 
are seen to increase with Ix/Iy, 12 /Iy, and R 2 max* region 

1 - Iz/Iy 0 for Ix 0, reorientation costs approach (iT/4)R2jav 

derived in appendix C, The same costs apply also to the region (1 - 1^/ly) -*• 1 
where I 2 "*• 0. These results are plotted in figure 7. 
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Although two-impulse coning solutions were not obtained for general reori- 
entations of asymmetric bodies, such solutions can reasonably be anticipated 
because the three input quantities to the maneuver (coordinate axes torques) 
should allow control of three output quantities (final Euler angles) regardless 
of body inertial configuration. Such solutions may result in smaller maneuver 
costs than the three-impulse asymmetric results of figure 7. 

Figure 8 was prepared to illustrate the effect of R 2 max average cost 
for axisymmetric bodies. Average cost is seen to increase with maneuver range, 
particularly for small maneuver ranges. A similar plot for asymmetric bodies 
is shown in figure 9. Although these results are for the set of bodies having 
Ijj/Iy = 1.2, they are typical of other Ix^Iy sets. 

Thus far, results have been presented for general reorientations. Similar 
results are available for spin-axis reorientations. The effect of inertial 
distribution on the cost of spin-axis reorientations is illustrated in fig- 
ure 10 for maneuver ranges of tt, 11 / 2 , n/4, and 0.00873 radians. Reorientation 
costs for axisymmetric bodies (Ix/Iy =1) appear to be relatively insensitive 
to the inertia configuration. Costs for asymmetric bodies tend to be smaller 
for elongated bodies. As would be expected, a comparison of figures 7 and 10 
shows considerably smaller costs for spin-axis reorientations than for general 
reorientations . 

The cost relationship with maneuver range for spin-axis maneuvers of roll- 
symmetric bodies is shown in figure 11. By definition, the final roll attitude 
for spin-axis reorientations is arbitrary. This not only eliminates the need 
for a follow-up spin maneuver, but greatly increases the number of candidate 
coning solutions; both of these factors act to reduce maneuver costs signifi- 
cantly and in proportion to the spin inertia as evidenced by a comparison 
of figures 8 and 1 1 . 

Spin-axis maneuver costs for asymmetric bodies of the set Ix^^y “ ^ ^ 
are shown in figure 12 as a function of maneuver range for a range of values 
of 1 - 12 /Iy. Again, spin-axis maneuver costs are relatively small, increase 
with maneuver range, and are less sensitive (than general reorientations) to 
changes in I 2 because final attitude about the spin axis is arbitrary. 

The average cost results presented in this paper for asymmetric bodies 
represent minimum cost reorientations involving a mixture of polhodes about 
the X- and z-axes. Although polhode axes for minimum maneuver cost are unpre- 
dictable for individual reorientations, a decided preference for polhodes to 
surround the axis of smallest inertia is observed. An even stronger tendency 
is for polhodes to be about z for small R 2 ; however, as R 2 becomes larger, 
the likelihood for polhode about x increases. Similarly, the relative cost 
of inward compared with outward nutations is not predictable for individual 
maneuvers. However, it is noticed that inward nutations (opposed mainly 
by Ix) are generally associated with lower coning costs for bodies with 
Ix < ly and also that outward nutations (opposed mainly by ly) are usually 
associated with lower costs for bodies with Ix ^ ly 

Maneuver costs have been presented in the form of average cost for speci- 
fied maneuver ranges. A measure of how costs are dispersed about these average 
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values is given by their standard deviation. Ratios of standard deviation to 
average cost were computed for all inertial distributions and for maneuver 
ranges of 0. 00873 > ir/2, and ti radians. These ratios were found to be 

moderately constant within each class of inertial distribution and maneuver 
range. For this reason only the mean value r and standard deviation o of 
these ratios are presented in the following table as a function of maneuver 
range and type of body symmetry for both general and spin-axis reorientations: 


Maneuver 

Type of body symmetry 





range , 








"2 , max » 





radians 

Mean value 

Standard deviation 

Mean value 

Standard deviation 


of ratio, r 

of ratio, a 

of ratio, r 

of ratio, o 

General reorientations 

0.00873 

0.572 

0.004 

0.430 

0.024 

n/4 

.438 

.074 

.341 

.020 

T7/2 

.379 

.041 

.320 

.016 

TT 

.341 

.041 

.320 

. 

.032 

Spin-axis reorientations 

0.00873 

0.358 

0.002 

0.575 

0.023 

■jr/4 

.368 

.006 

.433 

.058 

TT/2 

.385 

.006 

.367 

.030 

IT 

.417 

.008 

.366 

.041 


CONCLUDING REMARKS 

This report presents data figures which can be used to determine average 
costs of reorientation by impulse coning for bodies of all possible inertial 
distribution and over a wide range of maneuver angles for spin-axis reorienta- 
tions as well as general reorientations. Impulse coning is a method of atti- 
tude reorientation in which impulsive torques are used to initiate and later 
to terminate free precessional motion of a body. In spin-axis reorientations, 
the final rotational attitude about the spin axis is not specified. 

All attitude reorientations of axisymmetric bodies can be accomplished 
with two-impulse coning. Moreover, for any given reorientation, minimal two- 
impulse coning costs are always less than three-impulse (coning plus spin 
maneuver) reorientation costs. 

All general reorientations of asymmetric bodies achieved by the impulse 
coning method and constrained to have equal initial and final angles between 
the spin axis and the angular momentum vector are found to require a spin cor- 
rection maneuver at the end of coning. As a result, average costs of reorien- 
tation are considerably larger than those for axisymmetric bodies. However, 
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it is believed that the aforementioned constraint can be eliminated from coning 
reorientations of asymmetric bodies with the result that spin correction maneu- 
vers would not be needed and average costs of reorientation could be reduced. 

Spin-axis attitude reorientations of bodies of any inertial distribution 
can be accomplished with two-impulse coning maneuvers. Average costs of spin- 
axis reorientations are considerably less than similar costs for general reori- 
entations, especially for foreshortened bodies. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 15, 1977 
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APPENDIX A 


PROGRAM ASYMCON 

The computer program ASYMCON, used to determine minimum costs for attitude 
reorientation of asymmetric bodies by the impulse coning method, is presented 
in this appendix. 


Symbols 

Pertinent symbols used in program ASYMCON and their definition or relation- 
ship to symbols used throughout the report are as follows: 


A 

ACOST 


indicates polhode about x-axis 


weighted average cost of maneuvers within given R 2 max* 

^(Q3. mini sin R 2 I) 

Z( I sin R 2 I ) 


C 

DEL 
FLAG 
IX,IY,IZ 
PI 
QR 

RF 

R1,R2,R3 

R1C,R2C,R3C maneuver about x polhode axis equivalent to R-|,R 2 ,R 3 maneu- 
ver about z polhode axis 

R12 Qq 

R12A 0av 

R12M extremum value of 0 

R18 F 

R20 A<j) 


indicates polhode about z-axis 
00 

FLAG 

I^, ly, and I 2 , respectively 


ir 

tg/tf 

0 . 5^1 + 


I'Pf^l 


1^2, A 


R-|, R 2 , and R 3 , respectively 
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R201 

latest value of A<|) 


R21 

‘t’o 


R26 

Q3,min 


R29 

e 


R31 

Q2 


R32 

2|e|lz/Iy 


R4,R5,R6 

Ix, ly, and I 2 , respectively 


R40.R41 

lower and upper sample limit for R-| , respectively 


R42,R43 

lower and upper sample limit for R 2 , respectively; 

R43 

R44.R45 

lower and upper sample limit for R 3 , respectively 


R54 

absolute value of sin R 2 ; used as weighting factor 
average cost 

for ^ 

R7 

'I'f 


R9 

'I's 


SD 

standard deviation of average cost 





SGN7 



SR4 

temporary storage for Ix or during polhode change 

TCOST 

weighted sum of maneuver costs 


TC2 

weighted sum of squared maneuver costs 


TNO 

weighted sum of maneuvers 
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Program Listing 

PROGRAM ASYMCON ( I NPUT ^OUTPUT » TAP£5= I NPUT t TAPE6=OUTPUT ) 

common R1 .R2.RA«R5»R6*R7.R9.Rl2tRia«R20.PI . FL AG « R 1 2A . SGN7 « R2 1 .R201 

COMMON CHANGE *R35 «R37 

real lX.IY.IZ 

real Ml 1 .MOl 1 .MFl 1 ,M12,M0l2tMF12«M13«M013.MF13 
real M21 .M021 ‘MFai ,M22.M022.MF22.M23» MF23 

real M31 .M031 •MF31 «M32 t M032 * MF32 . M33 « M033 . MF33 
integer CSAV 

integer A.C.Al ,CA1 

equivalence ( I X «R4 ) , ( I Y , R5 ) » ( I Z»R6 ) 

DA^a A.C/1HA*1HC/ 

NAMEL I ST/ I CONE/ 1 X « I Y, I Z « CON . PTDEN , R4 0 « R4 1 « R42 t R4 3 * R44 , R45 « PR I NT , SR 
PI =2.»ASIN< I . ) SR8= 1 80./P I SCON= 1 . S PRINT=1. * 1Y=100. S SR=0. 

PTDEN=8./PI * C0N=2. * R40=R42=0. S R44=-Pl S R4 1 =R43=R48=P I 
PTDEN=16./PI i R43=PI/2. 

1000 REAd(5. ICONE) 

IF (EOF ( 5 ) ) 1 4 . 1 5 

14 STOP 

15 CONTINUE 

R70=REF=0. S C26=1.E8 S TN0= TCOST = TC2=0 . S ISAV=0 * FLAG=0. 

KOUNT=0 S START=-1.E-6 S TCREP=0. 

R1 =pRi rH51 =R40+l ,/(2.*PT0EN) 

R2=PR2=R42+1 ./(2.*PT0EN) 

R3=PR3=R53=R44+CON/ <2.»PTDEN) 

WR I TE ( 5* I CONE ) 

R54=ABS(SIN(R2 ) ) t GO TO 2 

7 R3=PR3=R3+C0N/PT0EN S I F ( R3 . GT • R45 ) GO TO 25 S GO TO 2 

8 R1 =PR1 =R1 +1 ./PTDEN S 1 F { Rl . GT.R4 1 )GO TO 24 S GO TO 2 

9 R2=PR2 = R2+1 «/PTDEN S R54 = ABS ( S I N ( R2 ) ) S I F ( R2 . GT . R43 ) GO TO 26 

C compute EULER COSINE MATRIX FOR MANEUVER R1»R2*R3 ABOUT Z Y Z AXES 

2 SR1=SIN(R1) S SR2=SIN(R2) S SR3=SIN(R3) 

CRl=COS(Rl) S CR2=COS(R2) S CR3=C0S(P3) 

B1 1 =CR1 »CR2»CR3-SR1*SR3 S B1 2=SRl #CR2*CR3+CR1 *SR3 
B13 = -SR2*CR3 S 82 1 = -CR'l *CR2*SR3-SR 1 *CR3 
B22=-SR1*CR2*SR3+CR1*CR3 % B23=SR2*SR3 
B31=SR2»CR1 $ B32=SR2*SRl S B33=CR2 

C compute EULER ROTATIONS R1C« R2C*R3C ABOUT -X Y -X AXES 
Rl C=ATAN2 (B1 2«-Bl 3 ) S X3=-B33*S I N { R 1 C ) -B32*C0S ( R 1 C ) 
Y3=B22*C0S(R1C)+B23*SIN(R1C) S X2 =B2 1 *X3+B31 *Y 3 S Y2=B11 
R2C = ATAN2 (X2 »Y2 ) $ R3C = ATAN2 ( X3 ♦ Y3 ) 

3 IF (REP.EQ-0. )G0 TO I S R1=-R1C S R2=R2C i R3=-R3C 
C initialize FOR NEWTON RAPHSON ITERATION METHOD 

1 CHAnGE=0» S SGN7=0, S if ( I SAV.GT • 1 , ) SGN7= 1 . S K0UNT=0 
Q= (R4-R6 )*( .5-SGN7) $ R20=S61=-ABS(Q)/Q*START 
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CALL ASYITER S I F ( CHANGE . EQ. 1 • ) GO TO 18 
S71 =R20 1 

R20=S62= (R20+R20J )/2. 

CALL ASYITER $ IF (CHANGE.EQ. 1 • >GO TO 18 
S72=R20 1 

lOl R1 00= (S72-S71 )/(S62-S61 ) 

R20=S63= (S7I -Rl 00*S61 »/( 1 .-RIOO ) 

KOUNT=KOUNT+l S I F ( KOUNT , EQ, 1 1 ) GO TO 18 
CALL ASYITER * I F < CHANGE.EO. 1 • ) GO TO 18 
S73=R201 

C ACCEPTANCE CONDITION FOR ITERATION 

IF ( .001 /R8.GT.ABS ( (R20-R201 )/R20 ) )GO TO 31 
S61=S62 $ S62=S63 * S71=S72 $ S72=S73 
GO TO 101 

31 R21 =-Rl +PI/2.+R9/2.+PI *SGN7 

c compute costs of maneuver 

R22=ABS (2,*S IN (Rl2 ) *R7*R4/R 1 8 ) /R5 
R23 = ABS <2,*C0S (Rl 2 )*R7*R4/RI8 )/R5 

R29 = R1+R3-R9-R20 S R29 = R29-2 ,*P 1 * I NT ( R29/ ( 2 • *P I ) ) 
R31=R22*(ABS(SIN(R21 ) ) +ABS ( COS ( R2 1 ) ) +ABS ( S 1 N < R20+R2 1 ) )+ABS 
1 (C0S(R20+R21 ) ) 1/2.+R23 

C MINIMIZING COST OF SPIN CORRECTION MANEUVER 
R32 = 2.*-R6/R5*ABS (R29 ) 

RF = R23*.5*(1 .+ABS(R7/R29)/’(R7/R29) ) * OR = SQRT ( R32/ ( R3 1 -RF ) ) 

IF <QR.GT,R32/R23 )QR=R32/R23 S R26= ( R3 1 -RF ) * ( I , +QR ) +R32+R32/QR 
A29=R29-2.*PI*ABS(R29)/R29 S R32 =2 . *R6/R5* ABS ( A29 ) 

RF = R23» ( 1 .+ABS (R7/A29 )/ (R7/A29 ) ) S AQR=SORT (R32/ (R31 -RF ) ) 
IF(AQR.GT.R32/R23)AQR=R32/R23 i A26= (R31-RF >*( 1 . +AOR ) +R32+R32/A OR 
IF (A26.GE.R26 )G0 TO 60 S QR=AOR $ R26=A26 S R29=A29 
60 R26=R26* ( 1 ,-SR 1+R31 *SR 
C OOLHODE CONDITION DETERMINED 

A1=A S IF ( (REP-.5 )*R35/R37.LT.O, )A1 =C 

C lowest cost solution is SAVED 

IF (R26.GE.C26 )G0 TO 18 

C7=R7 S C9=R9 S C20=R20 S DEL=R12 S C12A=Ri2A S CAl=Al 
C26=R26 S C29=R29 $ C31=R31 * C2l=R21 i CSAV=ISAV 

CREP=REP » CCRl=Rl S CCR2=R2 S CCR3=R3 S CQR=QR 
CSTART=-ABS (0)/Q*START 

18 ISAV=ISAV+1 $ FLAG=0, $ IF ( I SAV.EQ. 1 >FLAG= 1 . 

I F ( ISAV.EQ.3 )FLAG=1 . $ I F ( I SAV.LT ,4 ) GO TO 1 S ISAV = 0 
C SETTING CONDITIONS FOR'POLHODE ABOUT ( A ) AXIS OR RESETTING FOR 
C POLHOOE ABOUT (C) AXIS 

6 SR4=R4 S R4=R6 S R6=SR4 

IF (REP.GT.O. )GO TO 5 S REP=REP+1 • S GO TO 3 
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5 REP=0. S Rl=P«l $ R2=PR2 S R3=PR3 
C resetting initial value of R20 FOR ITERATION 

IF (ABS ( START ) .NE. 1 .E-6 )GO TO 11 i START = -.‘5 S GO TO 3 
11 IF (C26 .lt .1 .E8 )GO TO 13 S START = -2 . *ST ART 

IF(A8S(START ),lT,5. )GO TO 3 S CSTART=-ABS (Q ) /Q*START S GO TO 16 
C COMPUTE STATISTICALLY AVERAGE MANEUVER COST 

13 tno=TNO+R54 S TC05T=TC0ST+C26»R54 S ACOST=TCOST/TnO 
TC2=TC2+C26*R54»C26 
TCREP=TCREP+CREP 

R70 = R70+1. S IF (PRINT. EQ.O. )GO TO 27 

c matrix transformation check for each maneuver 

10 S20=S2F=SIN<DEL ) i C20=C2F=C0S ( DEL ) S S30=SIN(C2l) 

C30=C0S(C21) S S1F=SIN(C7) S ClF=COS(C7) S S3F =S 1 N ( C20+C2 1 ) 

C3F = COS (C20 + C2 1 ) i Sl=SIN(CCRl) S C1=C0S(CCRI) S S2 = SIN(CCR2) 
C2=C0S<CCR2) S S3=SIN(CCR3~C29) £ C3=COS ( CCR3-C29 ) 

M011=C20*C30 S M03I =-S20*C30 

Ml 1 =-Sl »S3 + C 1 *C2*C3 £ M31=-S2*C3 

MFll =-SlF*S3F+ClF*C2F»C3F £ MF3 I =-S2F»C3F 
M012=-C20*S30 £ M032=S20»S30 

Ml 2=-Sl *C3-C 1 #C2*S3 £ M32=S2*S3 

MF12=-SIF»C3F-C1F*C2F»S3F £ MF32=S2F*S3F 

M013=S20 £ M033=C20 

M13=C1»S2 £ M33=C2 

MF13=C1F*S2F £ MF33=C2F 

M021 =S30 

M21 =C 1 »S3+S1 *C2*C3 . 

MF21 =C1 F*S3F+SIF»C2F*C3F 
M022=C30 

M22=C1 ♦C3-S1*C2*S3 
MF22=C 1 F*C3F-Sl F»C2F*S3F 
M23=S1 «-S2 
MF23=S1 F*S2F 

D1 l=M01 1*M11 +M01 2*M2 1 +M01 3*M31 -MF] I 
012=M01 1*M12+M012*M22+M0 i 3*M32-MF12 
D1 3=M01 1 *M1 3+MOl 2*M23+M01 3»M33-MF1 3 
02 1 =M02 1 *MI I +M022*M2 1 -MF2 1 
D22=M021 »M12+M022*M22-MF22 
D23=M021*M1 3+M022*M23-MF23 
D31 =M031 *M1 1 +M032*M2 1 +M033»M31 -MF3I 
D32=M03 1 *M1 2+M032 *m22+M033*M32-MF32 
D33=M031*M1 3+M032*M23+M033»M33-MF33 

IF(ABS(Ol 1 )+ABS(012»+ABS(013>+ABS<021 > +ABS ( 022 ) + ABS ( 023 ) + ABS ( 03 I ) 
1+ABS(D32)+ABS(033).GT..0001 ) 

1 WRITE (6. 53)01 1 .012.013.021 .022.023.031 .032.033 
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53 format { 1X9E1 1 .3/) 

16 CONTINUE 

WRITE(6«28)CCR1 . CCR2 . CCR3 » C7 tCREP ♦ CSA V * CST ART « C2 1 . DEL . C26 * C 1 2 A . 
lC9.C20,ACOSTtR70,C29,CAl ,PR1 .PR2.PR3.C0R 
28 FORMAT < 1 XF8.4 .F7. 4 .2F8.4 ,F3.0. I 1 ,F4.1 .F8.4,F7.4.F7.3.F7.4 *2F8.4 « 
1F5»2*F6»0»F8»4 < IXAl «F6«2«F5»2«F6»2*F5»2»F3»0 ) 

27 C26=1»E8 S START=-1.E~6 S GO TO 7 

25 R3=PR3=R53 $ GO TO 8 

24 R1=PR1=R51 S Sd=SQRT ( TC2/TNO-ACOST*ACOST ) $ SDA=SD/ACOST 
WR ITE ( 6 • 4 )R2 « TN0» TCREP « ACOST.SD*SDA 
GO TO 9 

4 F0RMAT(6X8E14.4 ) 

26 SD=SORT ( TC2/TN0-AC0ST*AC0ST ) $ SDA=ACOST/SD 

write <6*4 )R70«TN0»TCREP«AC0ST»SD«SDA 

GO TO 1000 
END 

SUBROUTINE ASYITER 

C THIS subroutine ITERATIVELY COMPUTES A SOLUTION SET OF VALUES FOR THE 
C VARIABLES INFLUENCING THE NUTATIONAL MOTION OF THT MANEUVERING BODY 

COMMON Rl,R2*R4«R5«R6*R7«R9«Rl2,R18»R20.PI.FLAGtR12At SGN7 *R2 1 .R201 
COMMON CHANGE «R35 «R37 

Ria= (R4+R5)/ (2.»R5)+ ( 1 • -2 . *FLAG > ♦ ( R5-R4 ) / ( 2 » *R5 ) *S 1 N ( R20 ) /R20 
R21 =-R20/'2. + PI/2.*FLAG 

R9 = 2.»(R2i+Rl'P!/2.~PI»SGN7) S R9 = R9-2 . *P 1 » I NT ( RP/ { 2 . *P I ) ) 

R12=ABS ( ATAN( TANIR2/2. )/SIN(R9/2. ) ) ) 

F1=TAN(R21) S F2=C0S(R21) * F3=TAN(R12) 

R35=R5* (R6-R4 ) +R4* ( R6-R5 ) *F 1 *F 1 

R36=R6» (R5-R4 )*Fl »F1 +R5* (R6-R4 )/ (F3*F3*F2*F2 ) 

R37=R6* (R4-R5 )+R4*(R6-R5 )/(F3*F3*F2*F2 ) 

IF (R35/R37.GT.0# )G0 TO 5 * CHANGE=1. S RETURN 

5 R12M=ATAN(SQRT(R35/R36) ) * 1 F ( FL AG . EQ . 1 . ) R 1 2M= AT AN ( SORT ( R35/R37 > ) 
RI 2A= (R12+R12M )/2. » 0 I FF = R I 2M-R 1 2 * D I FS=S I N ( R 1 2M ) -S I N ( R 1 2 ) 

IF ( ABS (DIFF ) .GT. ,001 ) R 1 2 A = ACOS ( 0 1 FS/D I FF ) 

R7=2»*ASIN(C0S(R2/2. )*SIN<R9/2. )/C0S(Rl2) ) 

C R7 PLACED INTO PROPER QUADRANT 

R50 = ABS (R9I/R9 i R5 1 = I NT ( R9/P I ) +R50 * R5 I = I NT ( R5 1 /2 . ) 

IF (R5I /2. ,EQ, INT (R51 /2, ) )G0 TO 1 S R7=2 , *R5 1 *P I -R7 S GO TO 2 

1 R7=2«*R51 *PI +R7 

C COMPUTE LATEST VALUE FOR R20 

2 R201 =R7*C0S(R12A )♦ (R4/R1 8-R6 )/R6 
RETURN S END 
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INERTIA CONSTRAINTS 

The following derivation generates boundaries for identification of all 
possible inertial distributions. 

The class of all rigid bodies is divided into two parts. Consider for 
the first part those bodies having 


I 


X 



(B1) 


where ly is intermediate inertia and Ij. is spin inertia. From the 
definitions 

Ix = J^[m(y^ + z2)j 

ly = r[m(x2 + z^)] 

I 2 = 2[m(x2 + y2^ 


it is apparent that inertia about any one of the three axes cannot exceed the 
sum of inertias about the other two axes. Thus, it follows that 


Iv + I. 


(B2) 


From condition (B1), the following relationships can be obtained: 


1 SO 

I., 


Ix 

— > 1 

ly 


(B3) 


From condition (B2) 


T - I ^ T 
-'■X -^z - -‘■y 


(B4) 
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Adding ly to both sides of inequality (B4) and dividing by ly give 


Ix / Iz , 

— + 1 U 2 


(B5) 


Equations (B3) and (B5) define triangular zone A of figure 6 where 
I3; S ly i 1 ^. 

For the second part, consider all remaining rigid bodies. These are 
restricted to those bodies having 


Ix ^ ly ^ Iz 


(B6) 


Applying the reasoning just preceding inequality (B2) to the second part yields 


Ix + ly ^ Iz 


(B7) 


Prom condition (B6), the following relationships can be obtained: 


1 




(B8) 


Dividing inequality (B7) by 


ly and rearranging give 


+ 



(B9) 


Inequalities (B8) and (B9) define triangular zone B in figure 6 for which 
Ix ^ ly = Iz- 


Zones A and B include all possible inertial distributions providing ly 
is the intermediate moment of inertia. 
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APPENDIX C 


AVERAGE REORIENTATION COSTS FOR ROD-SHAPED BODIES 

The derivation of the average reorientation costs for rod-shaped bodies 
is presented in this appendix. 

Consider first the configuration I^ "*■ 0 and 1 - Iz/Iy "*■ 0, a rod-shaped 
body along the x-axis. Because I^ "*■ 0, x is the polhode axis and the rod 
precesses about a momentum vector located in the y-z plane with a half-cone 
angle of 9 w ir/2. The average cost of reorientation can be written 

2 (Impulse x Time)niean 2 

Average cost = K = — K( Momentum x Time)njg^n 


= Y = 2K R2,mean 

where K is a factor to account for reaction thrusters fixed relative to 
the body principal axes and . w is the amplitude of the angular velocity 
vector. 

Assuming that the direction of the maneuver relative to the transverse 
thrusters is uniformly distributed about the z-axis, 


2 ij 

K = (Isin <l)ol + Icos <t>ol)niean " if Jq = ~ 


Values of R 2 mean average cost are listed as follows for selected values 

^2, max used in the body of the paper: 


®2 , max 

ir 

TI/2 

TT/4 

0.00873 

^2 , mean 

lT/2 

TT/3 

0.5480 

0.00617 

Average cost 

if 

8/3 

1.461 

0.0164 


Finally, for the configuration 1 - Iz/Iy I and I^ = ly, the body is a 
rod along the z-axis. Although the inertial notations have changed relative 
to the previous configuration, the physical situation is identical as are 
the reorientation costs. 
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(b) Angular momentum components. 

Figure 1.- Orientation of body and inertial 
coordinate axis systems for Euler rotation 
sequence about body axes z,y,z, 

respectively. 



z 



Figure 2.- Coning geometry showing momentum components. 
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Figure 3»- Symmetric solution pairs showing 
cone geometry and alternative paths for 
positive and negative coning rates. 
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Integrated solution trajectories for example maneuver, 
showing inward and outward nutation classes. 







Figure 7.- Effect of inertia distribution on average cost of reorientation. 




Figure 7.- Continued. 



Figure 7.- Continued. 




Figure 7.- Concluded. 



Figure 8,- Effect of maneuver range on average cost for reorientation over 
range of 1 - Iz/I values for axisymmetric bodies. 




Effect of maneuver range on average cost of reorientation for 
asymmetric bodies witn Ix^^v = ^•2- 






Figure 10.- Effect of inertia distribution on average cost of spin-axis reorientation 




Figure 10.- Continued. 





Figure 10.- Concluded. 
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